Peculiarity of the Coulombic Criticality ? 
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We indicate that Coulombic systems could correspond to Wilson effective Hamiltonians similar 
' to that of the ordinary (nonionic) fluids but with a negative (y3*-coefficient. In that case, solving the 

^\ I "exact" renormalization group equation in the local potential approximation, we show that close 

0^ ■ initial Hamiltonians may lead either to a first order transition or to an Ising-hke critical behavior, 

' the partition being formed by the tri-critical surface. Hence the theoretical wavering encountered 

f— I , in the literature concerning the nature of the Coulombic criticahty may not appear senseless. 

PACS numbers; 61.20.Qg,05.40.+j,05.70.Jk,64.60Fr 
In a number of experimental, theoretical and computer simulation studies, the problem of Coulombic criticality 
has been addressed For some ionic fluids mean field critical behavior is observed experimentally. Other systems 
demonstrate Ising-type behavior, while still other ionic fluids exhibit a crossover from classical behavior to Ising-like 
' as the critical point is approached. At the theoretical level, the confused state of the subject has been well described 
J> ■ and clarified by Fisher and Stell 
I/"} Coulombic criticality may be studied within the restricted primitive model (RPM: equal numbers 7V+ and iV_ 

0^ of positive and negative hard spheres of equal diameter, d, with charges, ±ze, immersed in a structureless solvent 
- of dielectric permittivity e; z/e = 1 in what follows). If one generally agrees on the existence of a liquid-gas- like 
transition at low concentration and low temperature in RPM, it is still not evident what kind of critical behavior the 
^ _ model actually possesses. Mean-field-like or Ising-like critical behavior, crossover from mean-field to Ising behavior, 

\ tricriticality or first order transition: all these conclusions have been expressed and discussed [Q. 
"Y"; ' Our aim, in this letter, is first to compute the Wilson effective Hamiltonian corresponding to RPM at its critical 
' point and then to apply the renormalization group (RG) techniques to study it. Let us summarize the main lines of 
^ . the calculation of the effective Hamiltonian: the details will be published elsewhere |Q . 

I ' The Hamiltonian of the RPM (in units of fcsT) is a sum of the repulsive (hard-core) part and of the coulombic 
\ part: 
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' here (j)hc (i^ij) are the hard-core interactions, {fi} denote coordinates of the N = -f N- particles, fij = — r}, and 

^ ' p = N/fl is the density, where ft is the volume of the system. The coulombic interactions in Eq.(^ are written in 

■ - - ■ terms of the charge density fiuctuation amplitudes, Pfc = "-^ ^ The latter are expressed in terms of the amplitudes 

of density fluctuations of positive (a), ni — Q^i '^f=i ^^P {~*^^}' ^^'^ negative (&), n'j. = '^f=i ^^P {"^^'^j' 

particles. In Eq.(0) i'{k) = {Attc^ /UbT^ /A^, and the primes over sums denote that terms with i = j in the first sum 

and with fc = in the second sum are excluded. 

Using the Hubbard-Schofield scheme we map the ionic fiuid Hamiltonian into a Wilson effective Hamiltonian 
similar to that obtained for ordinary (nonionic) fiuids Q which, discarding the derivatives of the field (local potential 
approximation), has the following form: 



H = df 



i(V^(r))Vl/(^(r)) 



(2) 



with the "potential function" : 
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where 6^ = 4tt (Stt)^/^ p*i/3/r* (p* = p(p is the reduced density and T* — kBTtfi/e"^ is the reduced temperature). As 
for the ordinary fluid the coefficients U2n may be expressed in terms of the cumulant averages, (p^ ■ ■ ■ Pj; ) > 

\ ^ / cR 

where averaging is performed over the reference system having only hard-core interactions. First we calculate these 
quantities for the lattice-gas model and find that U2n — (— f )"^^ which yield the usual Sine-Gordon Hamiltonian for 
the RPM 1). To perform evaluation of the ojf~lattice coefficients U2n we use a symmetry of the RPM with respect 
to the hard-core interactions, the definitions for the correlation functions, |^,^, and express the coefhcients in terms 
of the Fourier transforms (taken at zero wave-vectors) of the "cluster" functions of the reference hard-sphere system. 
In obvious notations §] these read: /i2(l,2) = 32(1,2) - 1, /i3(l,2,3) = 33(1, 2, 3) - 32(1, 2) - .92(1, 3) - 52(2, 3) -I- 2, 
etc., where (7n(l, ■ ■ ■ ,n) are n-particle correlation functions. Using the relation between the gn+i and gn M, 



d 

XP^jrP^'gn = Pp^ 
dp 



ngn + p dfn+l (gn+l - gn) 



(4) 



where x — p^^dp/dP is the compressibility, we iteratively express the Fourier transforms of hn at zero wave-vectors, 
/i„(0), in terms of /i„_i(0) and its density derivative, and ultimately in terms of /i2(0) and its density derivatives. Then 
we use the relation Q ph2{0) = pkBTx — 1 = -Zq, and obtain coefficients for the off-lattice effective Hamiltonian. In 
particular, U2 = 1, M4 = — (1 -I- 3zq), uq — 1 + 15(zq -|- zqZi + zi), . . . where zi = p (dzo/p) , ... To obtain zq one can 
use the virial expansion for the hard-sphere pressure, P/ pksT =1-1- V, Bkp'' with the coefficients Bi, . . . Bq known 
(for small densities), or the Carnahan-Starling equation of state |g,|6|. Therefore, applying the above scheme, all 
the coefficients of the effective potential V (ip) may in principle be found. 

We have studied the density dependence of U2n up to 2n = 14 and observed that all the coefficients are negative in 
the density interval ~ 0.07 ~ 0.09 where the critical density of the RPM is expected to be. We have then performed 
an empirical analysis and found that the boundaries of the density interval where the coefficients U2n are negative 
depend fairly linearly on 1/n (see Fig. 1). Extrapolating this dependence we have found that all the coefficients 
become positive for n > 22 (see Fig. 1). Being secure in the knowledge that the effective Hamiltonian is bounded 
from below, we can envisage a RG analysis. 

We thus consider the "exact" RG equation in the local potential approximation |^] in three dimensions which, using 
the same notations as in reads: 

in which y stands for the dimensionless field and f{y, I) — dV{y, l)/dy, f — df/dy, f" = d^f/dy^, / = df/dl with I 
the RG scale parameter (that relates two different "momentum" scales of reference such that A; — e^'Ao). 

Ideally, the question raised may be formulated as follows: considering the effective Hamiltonian for the RPM at its 
assumed critical point (taken e.g. from Monte Carlo data |^,|l^) as an initial Hamiltonian (1 = 0) for Eq. (^), will 
the solution of Eq. (||) flow toward the Ising fixed point [the unique non-trivial fixed point of Eq.(|^)] or not? 

Unfortunately, considering the function f(l = 0,y) which gives the initial conditions of RPM for Eq.(|^) (at p* = 
0.0857, T* = 0.052 §, or p* = 0.080, T* = 0.0488 0), we have observed that the denominator 1 + f'{y) in §) 
has singularities in the interval, 0.1 < y < 0.12. Moreover, Eq.(^) does not allow us to handle values of Hamiltonian 
coefficients as large as those of Eq.(^. Considering that Eq. (H) is an approximation, that impossibility has no 
particular significance relative to the behavior of RPM. 

However, we have considered that the negative value of U4 found for RPM could be a characteristic feature of 
qualitative importance for (some) ionic systems. Thus we turn our attention to the solutions of Eq.(^ with initial 
functions involving negative values of U4 . 

In j|] a detailed study of the approach to the Ising fixed point using Eq. (||) has been presented. However all the 
initial Hamiltonians considered were taken with U4 > 0. To our knowledge, the few studies based on RG techniques 
that have, up to now, considered negative values of U4, either perturbatively [pj]| in three dimensions or, more recently, 
non-perturbatively in four dimensions, have concluded that there is no stable fixed point (the Ising fixed point 
cannot be reached starting with U4 < 0) and thus to the lack of "true" criticality (there would be no divergent 
correlation length, like in a first order transition). 

To be short, we consider the following simple functions as initial conditions to Eq. (||) (a detailed study of the case 
W4 < will be published elsewhere p3[|): 

f{y, 0) = ^2(0)2; + ^4(0)2;^ -f ueiO)y' (6) 

corresponding to a point of coordinates (^2(0), U4(0), ue{0) ,0,0, • • •) in the space S of Hamiltonian coefficients (the 
dimension of S is infinite). Since we want to set U4(0) < 0, and at least one positive higher term is needed to make 
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the Hamiltonian bounded from below, we choose to set ue(Q) positive. Having chosen a (negative) value for U4(0) and 
a (positive) value for uq{0), we use the "shooting" method |^ to determine the critical value ^2(0) of ^2(0) which 
brings /(y,0) in the critical subspace Sc of S. The "shooting" method is based on the fact that, for sufficiently large 
values of I, the RG trajectories go away from 5c in two opposite directions according to the sign of ^2(0) — U2{0)- 
Let us summarize our results by considering the case 1*4(0) — —6 as an example (see Fig. 2). 

A If uq{0) = 16, we find 0.3836174 > ^2(0) > 0.3836151. The associated RG trajectory goes away from the Gaussian 
fixed point Pg and remains in the sector M4 < of 5c. Hence it never reaches the Ising fixed point that lies in the 
sector U4 > . Instead the trajectory is attracted to a stable submanifold of dimension one (an infra-red stable 
trajectory) that emerges from Pq. On a pure field theoretical point of view, this trajectory is a renormalized 
trajectory (let us denote it by r„^<o) to which is associated the continuum limit of an asymptotically free scalar 
field theory in three dimensions. It is very likelythat T„^<o is the continuation below four dimensions of the 
continuum limit recently studied on a lattice in |l2| and corresponding to a scalar theory with negative quartic 
interaction. The lack of any fixed point ending T„^<o means that the correlation length remains finite at the 
assumed critical point. This situation could be compared to the fact that no heat-capacity divergency has been 
observed in the Monte Carlo study of the RPM 

B If M6(0) = 20, we find 0.30131 > 1*2(0) > 0.30122 and the associated RG trajectory goes toward the (Wilson-Fisher) 
Ising fixed point and approaches it along the usual renormalized trajectory associated with the continuum limit 
of the scalar field theory in three dimension usually called the iy9|-field theory. This renormalized trajectory 
(denoted by Ti in ||]) interpolates between the Gaussian and the Ising fixed points. Hence there exist initial 
Hamiltonians with W4 < that belong, nevertheless, to the basin of attraction of the Ising-like fixed point. 

C Between the two preceding cases, we find a trajectory (with uq{Q) = 18.3125 • • • and 0.3324573 > 1*2(0) > 0.3324549) 
that directly flows towards Pq (it is neither attracted to T„_j<o nor to Ti). That kind of initial Hamiltonian 
obtained by adjusting two coefficients (1*2(0) and 1*6(0)) lies on the tri-critical subspace St of S. Any trajectory 
on St approaches Pg along a unique (attractive) trajectory that imposes the required very slow (logarithmic) 
flow in the vicinity of Pg (see Fig. 2). 

So, it appears that for 1*4(0) < 0, very close Hamiltonians may lead to very different behaviors and this feature 
is due to the vicinity of the tricritical subspace. If one adopts the idea that effective Hamiltonians for Coulombic 
systems may be characterized by a negative value of u^, then it is not amazing that the tricritical state has been 
the center of a recent discussion relative to the nature of the Coulombic criticality. Also, it is easy to verify that the 
experimental observations of mean field values for the critical indices might be due to a retarded crossover towards the 
Ising behavior . One may observe on Fig. 2 that the trajectory corresponding to case B follows a long path to the 
Ising fixed point (and passes close to Pg) compared to a trajectory that could correspond to a simple fiuid. In order 
to illustrate qualitatively that retarded crossover, we have drawn in Fig. 3 the evolutions [calculated from Eq. (|^)] of 
a pseudo effective exponent fog (r — |ii2(0) — ii2| /| 1*2!) associated to two effective Hamiltonians with 1*4(0) = —6 and 
M6(0) = 20 on one hand [RPM-like ? (B)] and 1x4(0) = +3 and W6(0) = on the other hand (Simple fluid). 

Helpful discussions with H. Stoof, G. Stell and B. Lee are gratefully acknowledged. The work was partly supported 
by Natural Science and Engineering Research Council of Canada. 
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FIGURE CAPTIONS 



Figure 1 Dependence of the boundaries of the density interval where the coefficients U2„ of the effective Hamiltonian 
are negative as a function of 1/n. Extrapolation suggests that all the coefficients with n > 22 are positive. 
(Note that the density at which the "negative" interval shrinks to zero, p — 0.0856 is very close to the critical 
density from the Monte Carlo data |^,0). 

Figure 2 Projection onto the plane {u2, U4} of various RG trajectories (in 5c) obtained by solving Eq. (||). Black cir- 
cles represent the Gaussian (Pq) and Ising (IFP) fixed points. The ideal trajectory (dot line) which interpolates 
between these two fixed points represents the renormalized trajectory (RT) of the so-called 0| field theory in 
three dimensions (usual RT). White circles represent the projections onto the plane of initial (unrenormalized) 
critical Hamiltonians. For ^4(0) > 0, some effective Hamiltonians run toward the Ising fixed point asymptot- 
ically along the usual RT (simple fluid). Instead, for U4(0) < and according to the values of Hamiltonian 
coefficients of higher order, the RG trajectories either (A) meet an endless RT emerging from Pq (dashed curve) 
and lying entirely in the sector U4 < or (B) meet the usual RT towards IFP. The frontier which separates 
these two very different cases (A and B) corresponds to initial Hamiltonians lying on the tri-critical subspace 
(white square C) that are source of RG trajectories flowing toward the Pq asymptotically along the tricritical 
RT. Notice that the coincidence of the initial point B with the RG trajectory starting at point A is not real (it 
is accidentally due to the projection onto a plane). The restricted primitive model could correspond either to 
the case A or to the case B. 

Figure 3 Sketchy representation of the effective exponents i'cff(T) calculated in the local potential approximation 
along the two RG trajectories which, on Fig. 2, both flow towards the Ising fixed point ["RPM-like ? (B)]" 
corresponds to ^4(0) < and "Simple fluid" to 1/4(0) > 0). The full lines roughly indicate the experimentally 
accessible parts in each case. 
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Signs of RPM Hamiltonian coefficients 
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